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Abstract: We have developed a realistic, fully general relativistic computer code to simulate optical 
projection in a strong, spherically symmetric gravitational field. The standard theoretical analysis of 
optical projection for an observer in the vicinity of a Schwarzschild black hole is extended to black hole 
spacetimes with a repulsive cosmological constant, i.e, Schwarzschild-de Sitter spacetimes. Influence 
of the cosmological constant is investigated for static observers and observers radially free-falling from 
the static radius. Simulations include effects of the gravitational lensing, multiple images, Doppler and 
gravitational frequency shift, as well as the intensity amplification. The code generates images of the sky 
for the static observer and a movie simulations of the changing sky for the radially free-falling observer. 
Techniques of parallel programming are applied to get a high performance and a fast run of the BHC 
simulation code. 
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1 Introduction 

General relativistic deflection of light and lensing effects in gravitational field of stars 
were firstly investigated by Einstein pE]. In the vicinity of relativistic compact objects 
(black holes or neutron stars) these effects have strong influence on properties of the 
optical projection which become different than those of the optics in the flat spacetime 
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as we experience it in our everyday life [3]. Several authors have developed ray-tracing 
or simulation computer codes for modeling general relativistic optical projection in the 
vicinity of rotating or non-rotating black holes and neutron stars without presence of a 
cosmological constant, see, e.g., [2l l5l l6l IHl l9l flOl [16] . 

Recent observations indicate the cosmic expansion accelerated by a dark energy that 
can be described by a repulsive cosmological constant, A > [HHT]. We investigate the 
influence of A > on the appearance of distant universe for observers in close vicinity of 
nonrotating Schwarzschild-de Sitter (SdS) black holes. In order to obtain a good qualita- 
tive picture of the A influence, our simulations have been performed with unrealistically 
high values of A. 



2 Schwarzschild— de Sitter geometry 

The line element of the SdS spacetime has in the standard Schwarzschild coordinates and 
geometric units (c = G = 1) the form 

ds 2 = _ ( 1 _ M£ _ AA d<2+ / 1 _2M_ - 1 ^ + r2 2 + gin2 Q d(j)2 
V r 3 J V r 3 / 

where M is mass of the central black hole, and A ~ 1CT 56 cm -2 is the repulsive cosmo- 
logical constant. It is advantageous to introduce a dimensionless cosmological parameter 
y by the relation y = |AM 2 . The location of horizons is given by the condition g u — 0. 
Two event horizons exist for y G (0,y CT i t ), where y CT \ t = 1/27. The black hole and the 
cosmological horizons are located at 

2 ?r + £ 2 tt-£ 
r h = —=cos , r c = —= cos , (2) 



respectively, where 

f = cos" 1 3\/%. (3) 

The spacetime is dynamic at r < rh and r > r c . The static radius, (a hypersurface 
where the gravitational attraction of the black hole is balanced by the cosmic repulsion) 
is located at 

r s = ?ri (4) 

With increasing value of y, the horizons aproach to each other. In the critical case of 
y = Z/crit = 1/27, the horizons and the static radius coincide at = 3. For of y > 1/27, 
the spacetime is dynamic at r > 0, and describes a naked singularity |13j . We consider 
only spacetimes admitting existence of static observers that have y < 1/27. 



3 Optical projection in Schwarzschild— de Sitter spacetimes 

Construction of relativistic optical projection consists of finding all null geodesies con- 
necting the source and the observer, i.e., solving the so-called emitter-observer problem. 
An observer will see the image generated by the concrete geodesic in direction tangent to 
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the photon trajectory in observer's local frame, therefore given by space part of locally 
measured 4-momentum of photons p^ s . Directional angle a related to the outward ra- 
dial direction and the frequency shift g of the photon (the ratio of observed and emitted 
energy) are given by the general relations 

M (*) 
cosa = -^f, (5) 

Pobs ^source 

The indeces 'obs' (observer) and 'source' denote the components locally measured by an 
observer or a source. In the SdS spacetimes null geodesies are characterised by the impact 
parameter b defined as the ratio of constants of motion, b = $/£ [J3J HI] . For an observer 
located at r Q b s , a and g are functions of b only, as shown in Appendix. 

Due to spherical symmetry of the SdS geometry (JTJ) , it is sufficient to consider only 
sources and observers located in the equatorial plane. Considering observers located at 
= 0, A0 along geodesies connecting the source and the observer reads 

A0 = -^ SO urce ~ ^TT, (6) 

where (j) source is an angular coordinate of the source. The image order k takes values of 
0, 1, 2, . . . , +oo for geodesies orbiting the central black hole clockwise, and —1, —2, . . . , — oo 
for geodesies orbiting the central black hole counter-clockwise. The first direct and indi- 
rect images correspond to k = and k = —1, respectively. 

Photon motion in the SdS spacetimes is governed by the Binet formula [U Q3] 

M ' 1 (7) 



du \Jb~ 2 — u 2 + 2u A + y 



where u = r~ l . The critical impact parameter, b c = y 27/ (1 — 27 y), corresponds to the 
circular photon geodesic, which is located at r p h = 3M for arbitrary value of A [TBI IT%] . 
Photons coming from distant universe with b < b c end up in the central singularity, while 
photons with b > b c return back towards the cosmological horizont [TJ H3] . Using the term 
under the square root in (j7]) as a motion reality condition, a straightforward calculation 
yields relation for a turning point r t of geodesies with b > b c : 



r t = = cos 

/3(y + 6- 2 ) 



- arccos ( — 3 y 3 (y + b 2 



(8) 



and a relation for the maximum impact parameter 6 max for an observer located at given 

fobs; 

&max = i ■ (9) 

\J U obs ~ 2 Kbs ~ V 

Therefore, A0 can be expressed as 

A0Kource, U ohs , b) = T / , — , (10) 

J u SOUIce y u — u -|- zu -)- y 
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for geodesies with b < b c , or for geodesies passing the observer position ahead of the 
turning point. For geodesies passing observer position beyond the turning point A(f> can 
be expressed as 

)— ± /"° bS ^ U =F /" tmn ^ U (11) 

In (II 01) and (TTIi) , the upper (lower) sign corresponds to geodesies orbiting clockwise 
(counter-clockwise). These integrals express A(f> along the photon path as a function 
F (b,u b s ,u source ,y). Equation (jUj) can then be rewriten in the following way : 

F (b, u obs , u source , y) + (p SOU rce + k2rr = . (12) 

This equation, that can be understood as an integral equation with an eigenvalue b de- 
termines b as an implicit function of the source and observer coordinates, the image order 
and A. Unfortunately, F (b,u b s ,u source ,y) can be expressed in terms of elliptic inte- 
grals only, and do not allow to obtain an explicit formula for b ((psource, k, u b s , u source , y). 
Consequently, the simulation code uses standard numerical integration and root finding 
methods. 



4 Numerical solution and code outputs 

Final equation ( [121) was numerically solved using the BHC code for optical projection 
of observers located in the static region between the horizons, as well as for observers 
located in the dynamic region under the black hole horizon. In the static region, static 
observers and observers radially free-falling from the static radius have been considered. 
In the dymamic region under the black hole horizon, where static observers cannot exist, 
the optical projection was constructed for observers radially free-falling only. We consider 
A = 5.10 -3 , A = 10 -2 and pure Schwarzschild case with A = 0. 

Left panels of Figs. [U El and [3] show b as a function of |A0| along the appropriate 
geodesic connecting distant source and observer located at r b s . The right panel of these 
figures show a for static and radially free-falling observers. If we consider spherical 
symmetry of the problem, equation © implies, that |A0| < 7r corresponds to the first 
direct images of distant sources, whereas larger |A</>| > n corresponds to images with 
higher order. 

As shown in Fig. [H for observers located above the circular photon orbit b increases 
up to a maximum value 6 max given by equation (jHJ), then decreases, and asymptoticaly 
aproaches b c from above. Values |A0| < \A(f>(b max )\ correspond to ingoing geodesies with 
b < b c , or to geodesies passing the observer position ahead the turning point. Values 
|A0| > \A(p(b max )\ correspond to geodesies passing the observer position beyond the 
turning point. Situation is different for observers located under the circular photon orbit 
where only geodesies with b < b c can exist. In this case b increases monotonously and 
asymptoticaly aproaches to b c from below. 
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In all the cases, a monotonously increases up to a maximum value, which determines 
the black region on the observer's sky. We determine the size of the black region as a 
function of r Q b s and A in the next section. 

Fig. H] shows samples of visualisation outputs of the BHC code, simulated optical 
projection of a well-known galaxy M104 "Sombrero" , virtually located behind the black 
hole on the optical axis. Images show typical lensing effects as the first and second 
Einstein rings, first direct image, first indirect inverted image and mergence of higher 
order images with the second Einstein ring around the black region. Another static 
images, as well as dynamic simulations, may be downloaded from our web site [12] . 




Fig. 1 Optical projection for observers located at the static region above the circular 
photon orbit at r Q b s = QM. Left panel: Impact parameter b as a function of A<f>. Right 
panels: Directional angles for static and radially free-falling observers. 
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Fig. 2 Optical projection for observers located in the static region under the circular 
photon orbit at r Q ^ s = 2AM. Left panel: Impact parameter b as a function of A<p. 
Right panel: Directional angles for static and radially free-falling observers. 
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Fig. 3 Optical projection for observers located under the black hole horizon at r Q b s = 
0.7M. Left panel: Impact parameter b as a function of A<p. Right panels: Directional 
angle for a radially free-falling observer. 




Fig. 4 Simulated appearance of M104 "Sombrero" located behind the black hole. Left 
panel: for a radially free-falling observer at r Q b s = 20M in a pure Schwarzschild case 
(with nondistorted image in the right-bottom corner). Right panel: for a static observer 
at r obs = 5M with A = 1(T 3 . 



5 Apparent angular size of the black hole 

The apparent angular size S of the black hole can be naturally defined as the observed 
angular size of the circular black region on the observer sky, in which no images of distant 
objects can exist, and only radiation originated under the circular photon orbit can be 
observed [H El [3]. For observers located above the circular photon orbit the boundary 
of the black region corresponds to outgoing geodesies with b approaching b c from above, 
while for observers located under the circular photon orbit the boundary corresponds to 
ingoing geodesies with b aproaching b c from below. 
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Fig. 5 Apparent angular size of the black hole as function of observer's radial coordinate. 
Left panel: static observers. Right panel: radially free-falling observers. 

In the case of static observers we have [lj 

S = 2arccosv4(r obs ,?/;6) , A(r, y; b) = ±]j 1 - ^ ^1 - ^ - yr 2 ^j . (13) 

Here + (— ) sign corresponds to observers located above (under) the circular photon 
orbit. Above the circular photon orbit increasing A causes downsizing of the black region, 
whereas under the circular photon orbit the black region grows with increasing A. In the 
limit case of observers located just on the circular photon orbit, S is independent of A. 
It is invariably tt, i.e., the black region always occupies just one half of the observer sky. 

In the case of observers radially free-falling from the static radius [U [EH E]? the 
apparent angular size of the black hole reads pQ 



Z(r obs , y) + - 3yWA(r oha , V, b)) 

S = 2 arccos ) (-, (14) 

y/l - 3yV^ + Z(r ohs , y)A(r ohs , y; b) 



where 



Z(r,y) = ^+yr 2 -3y 1 / 3 . (15) 

The A dependency is qualitatively different. For radially free-falling observers S grows 
with increasing cosmological constant at all values of the radial coordinate except the 
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central singularity, where S is invariably tt, similarly to the case of static observers located 
on the circular photon orbit. Consequently, the radially free-falling observer will always 
observe smaller S then the static observer at the same radial coordinate. 

6 Conclusions 

In this paper we discussed the influence of A > on the optical projection in strong, 
spherically symmetric gravitational field. The influence depends on the value of the di- 
mensionless cosmological parameter y. In the present universe with A ~ 10~ 56 cm -2 
values of y are y ~ 10~ 40 for stellar black holes and y ~ 10~ 25 for supermassive black 
holes in galactic nuclei. Observable effects can be expected for y > 10 -15 which corre- 
sponds to supergiant black holes with masses M > 10 15 M Q [15] . In the case of primordial 
black holes in the very early universe, with assumed high values of repulsive cosmological 
constant, one can expect even stronger effects. Considering the electroweak phase tran- 
sition at T ew ~ 100 GeV, we obtain an estimate of the primordial effective cosmological 
constant A ew ~ 0.028 cm -2 , while considering the quark confinement at T qc ~ 1 GeV we 
obtain A qc ~ 2.8 x 10 _10 cm~ 2 and consequently higher values of y [13] . 

BHC code generates numerical solutions of the governing equation of the projection 
and static as well dynamic visualization outputs. Results show peculiar influence of A 
on the apparent angular size of the black hole for observers in different local frames. 
This influence vanishes for static observers located at the circular photon orbit. For 
future studies we plan to extend our method and the BHC code in order to study axially 
symmetric spacetimes with repulsive cosmological constant. 
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Appendix: tetrads and directly measured quantities 

It follows from the central symmetry of the geometry ([T]) that the geodetical motion of 
test particles and photons is allowed in the central planes only. The existence of Killing 
vector fields £m and of the SdS spacetime implies the existence of two constants of 
motion 



and the photon motion is determined by the impact parameter 



(16) 




(17) 



The 4-momentum of the photon reads [H] 





(18) 
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where we introduce new variables 

B\r,y) = \- 2 --yr\ A{r, y, b) = ±Jl - B\r, y)^. (19) 

The + sign corresponds to photons receding from the black hole, while - sign corresponds 
to photons infalling into the black hole. 

In order to calculate directly measured quantities, one has to transform the 4-momentum 
of the photon into the local frame of the observer. The local components of 4-momentum 
for the observer at given r obs can be obtained using the appropriate tetrad of base 4- 
vectors , 1-forms co^ and transformation formulas 

uj( a ) = eJ^W , p (a) = e^V . (20) 



Static observers 

The static observers located at rest at r = const, 9 = const, = const are endowed by a 
local frame with an orthonormal tetrad of 1-forms [Hj 

u ® = B(r,y)dt, J r) = — r dr, u {e) =rd9, w w = rsin^d0. (21) 

B(r,y) 

The local components of 4-momentum of the photon moving in the equatorial plane are 
given by the relations [l z 



_(*) _ S _W _ A(r ohs ,y;b) w _ IS _ 

Pobs " B(r ohs , y) ' P ° bs ~ B(r ohs , y) ' P ° bs ~ r ~ r obs ' {ZZ) 

Using general formulas (jSJ), the directional angle and the frequency shift are given as [T4] 

m u\ B[r SOUIce ,y) /oo\ 

cosa stat = -A(r ohs , y, b) , g stat = —— — (23) 

B[r ohs ,y) 

Observers radially free-falling from the static radius 

Local components and tetrads for free-falling observers can be obtained using Lorentz 
boost between the local frames of the static observer and a moving one at given r Q ^ s . The 
orthonormal tetrad of 1-forms of appropriate local frame has the form 



uj® = Jl- 3y 1 / 3 dt + Z(r, y)B~ 2 (r, y) dr, (24) 



= Z(r,y)dt + ^Jl-3y 1 / 3 B- 2 (r,y)dr, (25) 
J §) = r d6, (26) 
= rsin6d(j), (27) 



where we introduced a new variable 



Z(r,y) = J~+yr 2 -3yi/*. (28) 
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The components of 4-momentum of the photon measured by a observer radially free- 
falling from the static radius at a given r b s are given by the relations [l l 



pl = B2 ^ Qhs y) ( V 1 - % 1/3 + ^obs, y)A(r ohs , y; &)) , (29) 
£ = B2{ f Qhs y) (z(r, y) + y/l - 3y^A(r ohs , y; &)) , (30) 

P± = ^ = ^ ( 31 ) 

'obs ' obs 

The directional angle and the frequency shift are given by the formulas [1 

(z(r ohs , y) + y/l - 3y^A(r ohs , y- b)) 

cosa fa ii = - ; = f , (32) 

/ l-3 y y s + Z{r ohs ,y)A{r ohs ,y;b) 



Pobs B(r source ,y) /qq\ 

fell = -7^ — - — , [oo) 

^source Z(r ohs , y) COS 01 + J 1 - 3I/ 1 / 3 
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